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ABSTRACT 

The class II masers of methanol are associated with the early stages of formation 
of high-mass stars. Modelling of these dense, dusty environments has demonstrated 
that pumping by infrared radiation can account for the observed masers. Collisions 
with other molecules in the ambient gas also play a significant role, but have not been 
well modelled in the past. Here wc examine the effects on the maser models of newly 
available collision rate coefficients for methanol. The new collision data does not alter 
which transitions become masers in the models, but does influence their brightness and 
the conditions under which they switch on and off. At gas temperatures above 100 K 
the effects are broadly consistent with a reduction in the overall collision cross-section. 
This means, for example, that a slightly higher gas density than identified previously 
can account for most of the observed masers in W3(OH). We have also examined the 
effects of including more excited state energy levels in the models, and find that these 
play a significant role only at dust temperatures above 300 K. An updated list of 
class II methanol maser candidates is presented. 
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1 INTRODUCTION 

Methanol, CH3OH, is an important constituent of the inter- 
stellar gas in star-forming regions. It has a rich microwave 
and millimetre spectrum as a result of its asymmetry, and 
low-energy torsional vibrations of the methyl group against 
the OH frame can be excited under interstellar conditions. 
Radioastronomical observations of methanol provide a use- 
ful probe of physical and chemical conditions. This is partic- 
ularly true in regions of high-mass star formation, where a 
large gas-phase abundance results from evaporation of solid 
state methanol from grain mantles. Furthermore, the phys- 
ical conditions in such environments give rise to masers, as 
well as thermal and quasi-thermal emission. The interpre- 
tation of all such emissions requires reliable models of the 
excitation processes. 

Until now, the biggest deficiency in models of methanol 
exciation has been the lack of reliable rate coefficients for 
collisions with other constituents of the gas. Most authors 
have either stuck to a Local Thermodynamic Equilibrium 
approach (assuming a Boltzmann population distribution), 
or have used propensity rules to approximate the collision 
rate coefficients, based on a small number of double reso- 



nance experiments by Lees & Haque (1974). Recently, the 
first accurately calculated rate coefficients for methanol exci- 
tation by collisions with He and para-H2 have become avail- 
able (Pottage, Flower & Davis 2001, 2002, 2004a, 2004b, 
hereafter PFD). These give the prospect of obtaining more 
accurate information from observations of methanol emis- 
sions in a variety of environments. 

Here we focus on the class II masers of methanol, which 
are a tracer of recent high-mass star formation. They are of- 
ten closely associated with OH masers as well as infrared 
and ucH n region continuum emission from such regions. 
They pinpoint the very early development of a high-mass 
star, and may in fact pre-date the development of other 
tracers. Their class I maser cousins are also found in star- 
formation regions, but are ususally well separated from the 
source of excitation. The most prevalent class II methanol 
maser line is the 5i — 60 A + transition at 6.668 GHz, known 
at some 500 sites (see catalogs of Malyshev & Sobolev 2003; 
Pestalozzi, Minier & Booth 2004). The 2 -3_i E transition 
at 12.172 GHz has also been detected at many of these loca- 
tions. There are around 20 further class II methanol maser 
transitions which have features in common with the 6.6- and 
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12.1-GHz masers; these are generally weaker and known in 
only a small number of sources. 

The class II masers can be accounted for by the model 
of Sobolev & Deguchi (1994), in which infrared radiation 
pumps methanol molecules to the second torsionally excited 
state. While the excitation is predominantly radiative, col- 
lisional quenching becomes significant at high densities. In 
this paper we revisit past work with this model (Sobolev, 
Cragg & Godfrey 1997a, 1997b; Cragg et al. 2001; Sutton 
et al. 2001; Cragg et al. 2004) to ascertain the impact of the 
new collision data. As well, we look at the effect of including 
more energy levels in the calculations, and specifically more 
torsionally and vibrationally excited states. 



2 MASER MODELLING 

In the spectrum of an interstellar molecule, the relative in- 
tensity in different lines reflects the physical conditions un- 
der which they were produced. At low interstellar densities 
the populations of the various energy levels usually do not 
follow a Boltzmann distribution, but are determined by a 
statistical equilibrium between excitation and de-excitation 
processes. Molecules move from one energy state to another 
by radiative (absorption, spontaneous or stimulated emis- 
sion of photons) or collisional means (due to interactions 
with other constituents of the gas, principally H2 and He). 
Methanol masers arise when the competition between these 
processes produces population inversions between the ir- 
regularly spaced energy levels, and when there is enough 
methanol present that amplification can develop along par- 
ticular lines of sight. In practice, when modelling these pro- 
cesses, the self-consistent calculation of level populations 
from the transition rates can only be done for a finite num- 
ber of energy levels, and if too few are inlcuded the pop- 
ulation figures become distorted. The models rely on the 
avalability of high-quality molecular data for the energies 
of the various accessible quantum states, for the Einstein 
coefficients or line strengths which determine the rates of 
radiative transitions, and for the collisional rate coefficients. 
The currently available molecular data for methanol is de- 
scribed in the next two subsections, followed by a description 
of the class II maser model. 



2.1 Energy levels of methanol 

Methanol molecules exist in two symmetry states, denoted 
A and E, which differ in the alignment of the nuclear spins. 
As the normal radiative and collisional processes do not 
interconvert these symmetry states, they are treated inde- 
pendently. The A and E species are usually assumed to be 
equally abundant, although an excess of the A species fol- 
lows from its lower ground-state energy if the molecules 
formed under very cold conditions (e.g. [A]/[E] = 1.44 if 
formed at 10 K). This is a possibility in star-forming re- 
gions where the bulk of the gas-phase methanol is believed 
to have evaporated from grain surfaces, where it formed or 
was deposited during an earlier evolutionary stage. 

The methanol eigenstates are labelled by the overall ro- 
tational angular momentum quantum number J. All radia- 
tive transitions are subject to the selection rule AJ=0,±1. 



Conventionally the states are also labelled by K, represent- 
ing the projection of J onto the axis of symmetry, although 
the asymmetry and the internal rotation mean that K is not 
a strict quantum number (as it would be for a rigid sym- 
metric rotor). The A symmetry species has close pairing of 
levels, which are labelled by a ± symmetry label (related 
to the parity quantum number) and nonnegative K, while 
the E species levels are instead labelled by a signed K (see 
Lees & Baker 1968). The strongest radiative transitions have 
AK=0,±1. Over 700 methanol rotational lines have been 
identified in interstellar surveys (Lovas 2004), with quan- 
tum numbers ranging from J=0 to J— 26. 

The lowest energy vibrational motion in methanol is the 
torsion or hindered internal rotation of the methyl group 
against the OH frame (t>t=l), which is excited at energy 
approximately 200 cm -1 or 300 K above the ground state. 
Over 200 lines from the first and second torsionally excited 
states (vt — 1,2) are known in warm interstellar sources such 
as Orion (Lovas 2004). This suggests either that the gas is 
warm enough that such high levels are populated by thermal 
means, or that molecules have absorbed infrared radiation 
at around 50 /im to populate the torsionally excited states. 
Radiative transitions between different torsional states with 
AK=±1 are strongly allowed due to the large degree of 
torsion-rotation interaction in methanol. The sequence of 
higher torsional states («t = l,2,3,4 etc) each has its asso- 
ciated manifold of rotational levels (see Fig0. The next 
vibrational mode is the CO stretch (wco=l), which appears 
at energy approximately 1000 cm -1 above the ground state, 
closely followed by the CH3 rocking mode. 

The infrared signatures of star-forming regions associ- 
ated with class II methanol masers typically indicate dust 
temperatures 100-200 K (De Buizer, Pina & Telesco 2000), 
and it seems likely that the gas temperatures are of the 
same order. Even without the infrared radiation, one may 
need energy levels as high as 400-600 K in the ground and 
first excited torsional states to calculate populations reli- 
ably at such temperatures. Too few energy levels will gen- 
erate truncation errors as the tail of the population distri- 
bution is artifically constrained. Masers can be very sen- 
sitive to such errors, as they stem from small differences 
between nearly equal populations. When infrared radiation 
from warm dust is included in the models, it acts to pro- 
mote methanol molecules to the torsionally excited states, 
from which they decay back to the ground state in one or 
more steps. This disrupts the pattern of population among 
the ground state levels, giving rise to inversions and hence 
the class II masers. Although the population which resides 
in the torsionally excited states is typically very small, the 
pumping and decay process plays a vital role in determining 
the populations of the low-lying levels. It is necessary to in- 
clude the same number of rotational levels in each accessible 
torsionally excited state so that all populated ground state 
levels may participate in the infrared pumping. Sobolev & 
Deguchi (1994) found that both the Vt=l and v t =2 levels 
were essential to the class II maser pumping. Here we inves- 
tigate for the first time the influence of even higher states. 

The largest set of energy levels used in our calculations 
is displayed in Fig. for E species methanol. Details are 
given in Appendix A. 
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Figure 1. Energy levels used in the largest maser model calculations for the E symmetry species of methanol. Left panel shows levels 
of the ground vibrational state, including rotational states up to an energy threshold with maximum J= 22, plus corresponding levels 
of the CO-stretch vibration. Right panel shows corresponding levels of the first 4 torsionally excited states. For clarity, vt = 2 and 4 are 
shown as heavier and narrower ladders compared with vt = 1 and 3. All levels are on the same scale and referred to the same ground 
state, with energy in units of cm -1 on left axis and K on right axis. 



2.2 Rate coefficients for collisions 

In the interstellar gas He atoms constitute about 20 percent 
of the collision partners, while the remainder is a mixture 
of ortho- and para-H2, with the para ground-state form pre- 
dominating at low interstellar temperatures. Lees & Haque 
(1974) undertook microwave double-resonance experiments 
to investigate collisions between E-species methanol with 
both He and H2 at room temperature. Their results led 
them to develop the following propensity rules, which have 
been used in excitation modelling to date in the absence of 
any other information. Collision rate coefficients favour the 
dipole selection rules (AJ=0,±1), and decrease as 1/| AJ| for 
|AJ|>1, while AK=0 transitions are favoured over AK=1 
transitions by a factor of 4. Our previous modelling of class II 
methanol masers employed these propensity rules in the for- 
mulation of Peng & Whiteoak (1993), extended to the A- 
species as described in Sobolev et al. (1997a), and applied 
to transitions within « t =0 and within Vt=l. 

Pottage et al. (2001) calculated rate coefficients for ro- 
tational excitation of methanol by He at temperatures 10 
and 20 K, for both A- and E-species methanol up to 3=7 
in the torsional ground state. Their results showed trends 
commensurate with the propensity rules of Lees & Haque. 
In subsequent work (Pottage et al. 2002) an improved treat- 
ment gave rate coefficients for temperatures between 5 and 
200 K, for methanol energy levels up to J=9 in the ground 
torsional state, plus data at temperature 20 K for the first 
excited torsional state. Cross-sections for torsionally inelas- 
tic collisions between v t =0 and vt=l were presented by Pot- 



tage et al. (2004a), again for collision partner He, but in this 
case the trends did not follow the propensity rules for tor- 
sionally elastic transitions. Rate coefficients for excitation 
by para-H2 were calculated by Pottage et al. (2004b), again 
for temperatures between 5 and 200 K, and for methanol 
energy levels up to J=9 but in the torsional ground state 
(vt=0) only. As yet there are no calculated rate coefficients 
for methanol interacting with ortho-H2, which remains a po- 
tentially significant source of uncertainty. 

Excitation modelling requires rate coefficients for col- 
lisional transitions between every pair of energy levels in- 
cluded in the calculation. Since methanol has a huge num- 
ber of levels which may be populated under interstellar con- 
ditions, even the pioneering work of PFD does not pro- 
vide enough data to fully solve the excitation problem. It 
does however provide data for the most important colli- 
sional transitions into and out of the energy levels involved 
in nearly all the observed class II masers (the exception be- 
ing the 23.1-GHz 92 — 10i A + maser). In our calculations the 
'missing' rates have been approximated by extrapolating the 
given rates in line with the propensity rules, as described in 
Appendix B. 



2.3 The Sobolev-Deguchi model 

In the model of Sobolev & Deguchi (1994) masers develop as 
the result of radiative pumping by infrared radiation from 
warm dust, described by dust temperature T d , filling factor 
Wd, and opacity Td(^/10 7 ) 2 , where v is the line frequency in 
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MHz. Radiative transfer is treated in the large velocity gra- 
dient (LVG) approximation, augmented by a beaming factor 
e _1 , defined as the ratio of the radial to tangential optical 
depths, to represent the elongation of the maser region along 
the line of sight (Castor 1970). The maser region is treated as 
a volume of gas of uniform density «h, kinetic temperature 
Tk, and methanol abundance Xm- The optical depth is pro- 
portional to the specific column density parameter Nm/AV, 
defined as the column density of A- or E-species methanol 
tangential to the line of sight, divided by the line width 
AV. In terms of these parameters the relative abundance of 
methanol becomes Xm = (2e _1 AV(Nm/ AV))/(nnL) when 
the A and E symmetry species are equally abundant, where 
L is the extent of the maser region along the line of sight. 
The model optionally includes free-free emission from an 
underlying ucHu region, which, if present, will be amplified 
by the masers. The H II background continuum has electron 
temperature T e and free- free turnover frequency / c , and is 
geometrically diluted by a factor Wini- Here we use Wd=0.5, 
r d =l, e _1 = 10, T c =10 4 K, and f e =12 GHz throughout, and 
explore ranges of the parameters Td, Tk, nu, Nm/AV and 
Wmi relevant to class II methanol masers. 



3 RESULTS 
3.1 Extra levels 

The number of levels required is a function of the model 
temperatures, since the gas temperature influences the pop- 
ulation distribution among the rotational levels, while the 
dust temperature mainly governs the distribution between 
the assorted torsional and vibrational states. The best way 
to demonstrate that enough levels have been included in a 
calculation is to repeat it using more levels and show that 
the differences are minor. 

Fig. displays model brightness temperatures for 4 se- 
lected class II methanol maser transitions: 5i — 6o A + at 
6.668 GHz, 2 - 3_i E at 12.178 GHz, 9 2 - 10i A+ at 
23.121 GHz, and 3i - 4 A+ at 107.013 GHz. The 6.6- 
and 12.1-GHz lines are the strongest observed masers, while 
the 107.0-GHz transition is the most common of the other, 
weaker masers. The 23.1-GHz transition involves higher en- 
ergy levels than any other of the observed class II masers, 
and also displays the most sensitivity to the number of lev- 
els included among the 23 maser transitions which we have 
examined in detail. The calculations displayed make use of 
the new PFD collision data, but the same trends are seen 
in calculations using the old propensity rule collision model, 
so the sensitivity of the 23.1-GHz maser is not due to the 
extrapolation of the PFD collision rate coefficients above 
J=9. The plots show the variation with gas temperature Tk, 
dust temperature Td, gas density iih, and specific column 
density of methanol Nm/AV. In each case the unvarying 
parameters are T k =150 K, T d = 175 K, n H = 10 7 cm -3 , 
and Nm/AV = 10 12 2 cm~ 3 s. These models do not in- 
clude ucH II background continuum radiation. Fig.|5|displays 
traces for 3 different energy level sets: the smallest has ro- 
tational levels up to J=18 for each of the torsional states 
v t — 0,1,2 while the largest has rotational levels up to J=22 
for u t =0,l,2,3,4 plus the vco~l vibrationally excited state. 
The intermediate set has rotational levels up to J=18 for 
vt=0,l,2,3. 



It is apparent from Fig. [5] that the effects of expand- 
ing the energy level set are generally modest. The smallest 
set produces masers which follow the same trends as those 
calculated from the largest set. The differences are most pro- 
nounced at gas temperatures Tk > 200 K, at dust temper- 
atures Td > 300 K, and at large methanol column densities 
Nm/AV > 10 13 cm _3 s. A detailed examination including 
further sets of levels shows that expanding the set of rota- 
tional energy levels from J=18 to J=22 generally has a lesser 
effect than expanding the set of torsional states from v t =2 
to v t =3. The intermediate set gives a close approximation 
to the largest data set except at high dust temperatures. 
At T d > 300 K the CO-stretch vibration starts to have a 
significant quenching effect on the masers. But as we have 
not included (for lack of available data) the CH3 rock vi- 
bration which is comparable in energy, we cannot be sure 
that this regime is reliably modelled even with the largest 
set. In the remainder of this paper we use the intermediate 
data set which includes rotational levels up to J=18 and 
torsional states up to ut=3, and we do not investigate dust 
temperatures above 300 K. 

3.2 New collisions 

Fig. |31 displays model brightness temperatures for 4 class II 
methanol maser transitions, as 4 of the model parameters are 
varied in the way described in the previous section. Dotted 
traces represent the results using the old propensity rule col- 
lisions, while the solid traces show the results with the new 
PFD calculated collision rate coefficients. Comparing Fig. 
with Fig. it is apparent that the new collision data has 
much greater impact on the modelling than including more 
energy levels. Nevertheless, the behaviour of the masers is 
recognisably similar with the old and new collision models. 
All 4 transitions displayed become masers of similar peak 
intensity with both collision models, but there are some- 
times significant changes in the conditions under which they 
switch on and off. The same is true of 19 other observed or 
potential maser lines which we have examined in detail. 

At low gas density (10 4 cm -3 ) the collision model is 
irrelevant and the masers are entirely governed by radiative 
effects. At high gas density (10 9 cm -3 ) all the masers are 
thermally quenched. The collisions are significant over the 
intermediate density range, as the masers become progres- 
sively quenched. With the new collisions at Tk=150 K, all 
masers extend to higher density than with the old. The den- 
sity required for thermalization is increased in some cases 
by an order of magnitude. Although some masers are af- 
fected more than others, the effects are consistent with an 
overall reduction in the collision cross-section with the new 
model. At Tk=150 K the masers appear at a lower dust 
temperature threshold, and extend over a wider range of 
specific column density, with the new collision model. These 
effects are seen when the gas density is moderately high 
(here nH = 10 7 cm -3 ), meaning that collisions contribute 
significantly to the balance of effects in the maser pumping. 
Both these effects are consistent with an overall reduction 
in collision cross-section, so that the radiative pumping be- 
comes more influential. 

The picture changes as the gas temperature varies. At 
gas temperatures above 100 K the new collisions produce en- 
hanced maser intensity at high density, as described above. 
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Figure 2. Effects of parameter variation on brightness temperatures for 4 masers, using 3 different data sets. Solid lines include J=0- 
22 for «t=0-4 and ^co = li dashed lines include J=0-18 with i>t=0-3, dotted lines include J=0-18 with i>t=0-2. Fixed parameters are 
T k =150 K, T d = 175 K, nn = 10 7 cm -3 , Nm/AV = 10 12 ' 2 cm _3 s, with no Hn region background radiation, and using the new PFD 
collision data. 



At gas temperatures below 100 K it is not so clear-cut. For 
example, at 21=30 K some masers extend to higher density 
with the new collisions, while others are thermalized at lower 
density than before. Thus the significant collision processes 
(whatever these may be) have become more probable for 
some transitions and less probable for others. Correspond- 
ing effects appear in the behaviour as specific column density 
and dust temperature are varied for Tk=30 K. The net effect 
is that at low gas temperatures some of the weaker maser 
lines are enhanced while others are diminished. 



3.3 Limits on most common maser sources 

The majority of class II methanol maser sources have been 
observed only at 6.6 GHz. Surveys at 107.0 and 23.1 GHz 
have demonstrated that the weaker maser transitions are 
rarely detected at 6.6-GHz maser sites (Val'tts et al. 1995, 
1999; Caswell et al. 2000; Minier & Booth 2002; Cragg et al. 
2004). It is evident from Fig.|3]that the model predicts emis- 
sion from the weaker masers simultaneously with the 6.6- 
GHz maser under many combinations of parameters. How- 
ever, the 6.6-GHz maser is the brightest and the most readily 
excited under the conditions investigated. Sources in which 
only the 6.6-GHz maser is detected may experience condi- 
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Figure 3. Effects of parameter variation on brightness temperatures for 4 masers, using different collision rate coefficients. Solid lines 
use new PFD collision data, dotted lines use propensity rules. Fixed parameters arc T k =150 K, T d = 175 K, njj = 10 7 cm -3 , Nm/AV = 
10 12 ' 2 cm~ 3 s, with no Ho region background radiation, and using energy levels up to J=18 and vt=3. 



tions near the edges of the maser regime, or may have maser 
emission in the weaker transitions which is below the de- 
tectable limit. In Cragg et al. (2004) we used the upper lim- 
its on maser emission at 107.0 and 23.1 GHz to investigate 
conditions in 28 class II maser sources detected at 6.6 GHz, 
by using these constraints to define models which give ratios 
compatible with the observations. We reported model cal- 
culations for gas temperatures 75, 125, 175 K, dust 
temperatures Td=125, 175, 225, 275 K, gas densities in the 
range 10 4 < nu < 10 9 cm~ 3 , and methanol specific column 
densities in the range 10 10 < N M /AV < 10 14 cm~ 3 s. The 
background ucH II radiation was omitted in this case, since 
less than half the class II methanol maser sources have as- 



sociated ultra-compact H II regions, and it is suggested that 
the majority of masers accompany the very early stages of 
high-mass star formation, before the development of a de- 
tectable ucHu continuum (Ellingsen et al. 1996). These cal- 
culations used the propensity rules to define collisional exci- 
tation rates, and we have now repeated them using the new 
PFD calculated collision rate coefficients. While the new re- 
sults may be considered more reliable, note that the extrap- 
olation of the PFD data for J>9 leaves some uncertainty 
over the behaviour of the 23.1-GHz 9 2 - 10i A+ maser. 

We find that as a general rule the 6.6-GHz maser ex- 
tends to greater values of gas density with the new collisions 
than with the old, sometimes as high as 10 9 cm~ 3 , consis- 
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Figure 4. Heavy and dashed contours delimit the model parame- 
ter combinations consistent with the strongest and weakest upper 
limits on maser emission from survey at 23.1 GHz, while light and 
dotted contours do the same for 107.0 GHz. The shaded ares rep- 
resent conditions satisfying upper limits for both transitions, for 
sources with masers detected at 6.6 GHz, but not at 23.1 and 
107.0 GHz. The line log(iV M /AV) < log(ra H ) + 6.5 distinguishes 
plausible model conditions (lower right) from those requiring un- 
realistically large maser path lengths or methanol abundance (up- 
per left). Top panel is from Cragg et al. (2004) using the propen- 
sity rules for collisions, while bottom panel is based on the PFD 
collision rate coefficients. 



tent with the examples in Fig. |31 The difference is greatest 
when 21=175 K and Td = 125 K, where the 6.6-GHz masers 
extend to more than an order of magnitude greater den- 
sity than previously. This suggests that relevant collisional 
processes which act to quench the masers have become less 
probable with the new rate coefficents. The same trend can 
be seen for the 23.1-GHz masers, whereas for the 107.0-GHz 
masers the effect is less pronounced at high gas temperatures 
and sometimes reversed at low gas temperatures. Consis- 
tent with this pattern, the masers often extend to slightly 
higher values of /Vm/AV and switch on at lower thresh- 
olds of Td than in previous models. The analysis of Cragg 
et al. (2004) was based on ratios of brightness temperature 
in the various transitions, since absolute values cannot be 
determined from observed flux densities without knowledge 
of the source size. Repeating this analysis to define param- 
eter combinations consistent with observed limits on maser 
emission, we obtain similar results to before, but generally 



shifted to slightly higher values of gas density, as the follow- 
ing example illustrates. 

Fig. 0] compares the range of uh and Nm / A V com- 
patible with the observed upper limits at 21=125 K and 
Td=125 K with the old and new collision models. Dark 
and light shaded regions show parameter ranges in which 
brightness temperatures of the 23.1- and 107.0-GHz masers 
relative to the 6.6-GHz maser are below the strongest and 
weakest limits (respectively) set by the survey data. The re- 
sults with the new collisions can be summarized as a shift 
to higher gas density by about a factor of 10, so that some- 
what higher gas densities are now allowed by the survey lim- 
its, together with some relaxation of the limit on methanol 
column density. A shift in the same direction is seen with 
other 21 and Td combinations. The change is particularly 
apparent when 21>21 (as in Fig.^J, because the restriction 
log(iVjy[/AV) < log(nn)+6.5 (which ensures reasonable val- 
ues of the methanol abundance Xm and maser path length 
L) now leaves a more substantial range of plausible models. 

3.4 W3(OH) 

There are a few exceptional class II methanol maser sources 
in which masers have been detected at many frequencies. 
Although they have complex spectra at 6.6 and 12.1 GHz, 
with many maser components spaced over a narrow range 
of velocities, usually only one or two components are evi- 
dent in the weaker maser lines. When the velocities match 
it suggests that the masers at different frequencies are simul- 
taneously excited in the same volume of gas. Interferomet- 
ric observations have confirmed that the 6.6- and 12.1- GHz 
emission in some sources is spatially coincident to within a 
few milliarcseconds (Menten et al. 1992; Norris et al. 1993; 
Minier, Booth & Conway 2000). In the models such simul- 
taneous excitation is readily achieved (in fact it is difficult 
to find conditions which would explain emission in one tran- 
sition alone), but the combinations of masers which appear 
together are sensitive to the conditions. Thus, when masers 
are detected in several transitions, they provide a sensitive 
probe of the local conditions, so long as they coincide. Upper 
limits on maser nondetections can also provide constraints 
on the physical conditions, and with fewer assumptions since 
the limits apply to all maser sites within the beam. 

In Sutton et al. (2001) we exploited this method to esti- 
mate physical parameters for the maser region in W3(OH), 
based on BIMA observations of maser spikes at 86.6, 86.9 
and 107.0 GHz at velocity -43.1 kms -1 , which were found 
to be located near the northern edge of the ucH II region. No 
corresponding spike emission was obtained for the 80.9, 84.4, 
85.5, 94.5, 108.8, or 111.2-GHz maser candidates. While 
maser emission at several other methanol frequencies has 
been detected at the same velocity (Tabled, there are at 
least 2 distinct maser regions in W3(OH) emitting at this 
velocity (Menten et al. 1992), and so it is not readily ap- 
parent from single-dish observations whether the remaining 
masers coincide. Indeed, the 23.1-GHz maser peak emission 
is from the southern edge of the ucHu (Menten et al. 1988), 
unlike the 107.0-GHz maser (Sutton et al. 2001), suggesting 
that local conditions at these two sites are rather different. 

We have undertaken further model calculations based 
on the new PFD collision rate coefficients to compare with 
the W3(OH) observations. These models include continuum 
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Table 1. Comparison between maser flux densities observed towards W3(OH) and 4 model calculations, following Sutton et al. (2001). 
Observational references (a) Sutton et al. (2001), (b) Menten et al. (1988) , (c) Wilson et al. (1993) , (d) Haschick, Baan & Menten 
(1989) , (e) Slysh, Kalenskii & Val'tts (1995) , (f) Menten et al. (1992) , (g) Moscadelli et al. (1999) , (h) Wilson et al. (1985) . 
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radiation from an underlying ucHu region with VKhii=0.2. 
The results are summarized in Tabled which includes data 
(maser detections or upper limits) for the 9 transitions ob- 
served by Sutton et al. and included in the fit, 10 other tran- 
sitions with comparable maser profiles (but not necessarily 
emanating from the same site), and 4 further masers with 
complex multi-component profiles. The observed flux densi- 
ties are compared with the predictions of 4 models. For each 
model, calculated brightness temperatures are converted to 
flux densities based on a common source size which matches 
the observed flux density at 107.0 GHz. 

Model B reproduces the calculations of Sutton et al. 
(2001) using the propensity rule collisions, with param- 
eters T k =150, T d =175, log(n H )=7.0, log(JV M /AV)=12.0. 
This gives a good fit to the 86-GHz masers, and meets the 
upper limits at 80.9, 84.4, 85.5, 94.5, 108.8 and 111.2 GHz. 
Model C has the same parameters but uses the new colli- 
sions. The 86-GHz masers are far too strong, and the lim- 
its at 80.9, 84.4 and 108.8 GHz are not satisfied. Model D 



also uses the new collisions, but has a higher gas density, 
log(riH)=7.5. It gives a good fit to the 86-GHz intensities, 
but still fails the upper limits at 84.4 and 108.8 GHz. This 
is consistent with the conclusions of the last section, where 
(to a first approximation, and at high gas temperatures) 
the effect of the new collisions is comparable to a den- 
sity scaling, such that higher densities are now required. 
Model E again uses the new collisions, with lower temper- 
atures T k =T d =145, and log(n H )=7.0, log(JV M /AV)=12.0. 
The equal dust and gas temperatures generate masers of 
lower brightness temperature than the other models shown, 
and so a correspondingly larger source size is required to 
match the observed flux densities. It gives a good fit to the 
86-GHz intensities, and meets all the upper limits except 
108.8 GHz, where it is an improvement on Models C and 
D. Note that this model makes predictions more consistent 
than the others with the 157-GHz Jo-J-i E maser observa- 
tions. These observations are not included in the fit, because 
it is not known whether they emanate from the same north- 
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ern source, but in either case the other models overestimate 
them by at least an order of magnitude. Thus a compara- 
ble, or perhaps better, account of the observations can be 
obtained with the new collisions, under conditions which are 
not vastly different. 

3.5 Maser Candidates 

Sobolev et al. (1997b, hereafter SCG97b) examined the 
model predictions under a range of circumstances, and pre- 
sented a list of class II methanol maser candidates. These 
included all the class II maser transitions known at the time, 
a few which have been identified in subsequent observations 
(e.g. 85.568, 86.615 and 86.902 GHz by Cragg et al. 2001; 
108.893 GHz by Val'tts et al. 1999), some lines which have 
been observed in maser sources but which have not been 
found to display maser characteristics, and many further 
lines which have yet to be investigated. Table [5] presents 
an updated list, based on calculations with the new colli- 
sion rate coefficients of PFD. The range of model conditions 
chosen differs somewhat from those presented in SCG97b 
following our experience in fitting models to observations. 

Table [5] tabulates the masers for 8 selected models. All 
models have dust temperature Td > 100 K and specific col- 
umn density of methanol 11.5 < A?m/AV < 12.0 in the 
range which generate significant class II methanol maser 
intensities. They include a range of gas temperatures and 
densities to illustrate the effects of the collision model. Al- 
though the majority of 6.6-GHz maser sites have no asso- 
ciated ucHu continuum, there is evidence that the weaker 
maser transitions are more readily detectable in sources in 
which the ucHu region has already developed (Ellingsen 
et al. 2004), and so all but one of the selected models in- 
clude background continuum radiation. The list contains all 
transitions which attain brightness temperature above 10 4 K 
in any of the selected models. Some masers require rather 
specific conditions to become active, while others are bright 
across the full range of models illustrated. To make these 
trends more apparent, brightness temperatures above 10 6 K 
are shown in bold, while blank entries mean T b < 10 3 K. The 
6.6- and 12.1-GHz masers are usually the strongest, both in 
observations and the models, and we refer to the other tran- 
sitions as the weaker class II masers. The detectability of 
these weaker masers is governed by their (unknown) size, 
which determines the flux density; e.g. Slysh et al. (1999) 
obtained a lower limit of 5 x 10 K for the brightness temper- 
ature of the strongest 107.0-GHz maser spots in W3(OH). 

Model 1 is an example with low gas density where the 
masers are influenced by radiative effects alone. Since the 
gas temperature is irrelevant in these circumstances, it is 
equivalent to Model 7 of SCG97b. The same masers are 
present in the list as previously, although their brightness 
is reduced sightly, which may be attributed to the inclusion 
of more energy levels in the current modelling (see Fig. 
Model 2 represents warm gas (Tk=150 K) of moderately 
high density (nrr = 10 7 cm -3 ), where the collisions have 
begun to quench out the masers (e.g. there are only 3 tran- 
sitions with Tb > 10 6 K, compared with 13 transitions in 
Model 1). Model 3 illustrates the cool gas (Tk=30 K), warm 
dust (Td = 175 K) combination which generates the largest 
number of intense masers (24 have Tb > 10 6 K). These con- 
ditions are identical with Model 6 of SCG97b. A detailed 



comparison shows that the same masers are present in both 
cases, some being a little stronger with the improved mod- 
elling while others are a little weaker, due primarily to the 
effects of the new collisions. The maser candidates listed 
for Models 1-3 are expected to be most readily detected in 
sources with masers at 107.0 GHz, which is the most preva- 
lent of the weaker masers. 

Models 4 and 5 illustrate conditions which may be more 
typical of class II maser sources, since they meet the ob- 
served upper limits on 23.1- and 107.0-GHz masers from 
Cragg et al. (2004) (see Section 3.3). Model 4 has no back- 
ground ucH n continuum while Model 5 includes it; the effect 
is to enhance the maser brightness temperatures, particu- 
larly at the lower end of the frequency range. The list of 
strong masers with Tb > 10 6 K is confined to the widely 
distributed 6.6- and 12.1-GHz transitions, plus the (difficult 
to observe) 57.032- and 60.531-GHz transitions. 

Model 6 is obtained from fitting millimetre wavelength 
maser observations in W3(OH), and is the same as Model E 
from Table Again it illustrates warm gas at moderately 
high density, and since in this case the gas and dust tem- 
peratures are equal, most of the masers present do not de- 
velop very strong brightness temperatures. Model 7 illus- 
trates a combination of masers similar to those observed in 
NGC 6334F, including intense masers at 19.9 and 23.1 GHz, 
generated by a warm dust/cool gas combination. 

Finally, Model 8 has dust temperature closer to the 
threshold for class II maser activity, and exhibits weak maser 
action in several transitions absent from the other mod- 
els. These include the 25-GHz J2-J1 E class I maser series, 
and transitions at 9.936, 20.171, 23.444 and 104.300 GHz 
which develop prominently in models with temperatures in 
between the class II and class I methanol maser domains. 
The transitions probably characterize a new 'intermediate' 
class of methanol masers, perhaps exemplified in W33-Met 
(Voronkov et al. 2004). Of these lines tracing intermediate- 
class masers, only the 9.9-GHz line has been searched for, 
and only in a very limited number of sources. Slysh, Kalen- 
skii & Val'tts (1993) detected a weak 9.9-GHz maser in W33- 
Met, but obtained no detections from a search of 5 class II 
maser sources. 



4 DISCUSSION 

A good example of the distortions which can arise from 
model calculations including too few energy levels is pre- 
sented by Sobolev & Deguchi (1994). In their Model Dl, 
which includes levels up to J=12 for torsionally excited 
states « t =0,l,2, the 12.1-GHz 2o — 3_i E transition be- 
comes a strong maser with Tb = 4.7 x 10 10 K under warm 
dust/cool gas conditions. However, Model Dla, which em- 
ploys the same physical parameters but omits the vt=2 lev- 
els, instead produces absorption in the 12.1-GHz line, with 
Tb = —1.1 x 10 4 K. Note that there is considerable overlap in 
energy between the vt = l and v t =2 levels (Fig.0. Thus the 
switch-on of the class II masers at dust temperatures above 
100 K is clearly attributable to excitation to the levels of the 
second torsionally excited state. In contrast, we have shown 
that further extending the energy level set to include levels 
of v t —3, «t=4 and the CO-stretch vibration (i>co=l) has a 
minor effect by comparison. We find the class II masers to 
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Table 2. Class II methanol maser candidate lines, which have brightness temperature > 10 4 K in one or more of the 8 models presented. 
All model calculations employ the new collision data of PFD. Blank entries signify brightness temperatures < 10 3 K, while brightness 
temperatures > 10 6 K are in bold. Entries which meet the specified thresholds, but represent quasi-thermal rather than maser emission, 
are marked with an asterisk. 
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Table 2 - continued 
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be adequately modelled for dust temperatures up to 300 K 
when levels up to J=18 and w t =3 are included. At higher 
dust temperatures excitation to the vco = l levels tends to 
quench the masers, and it is likely that the levels of the 
CH3 rock vibration will also play some role. The indications 
from infrared observations are that dust temperatures in 
class II methanol maser souces are in the range 100-200 K 
(De Buizer et al. 2000), consistent with these findings. 

The accurately calculated collision rate coefficients of 
PFD provide a much firmer footing for excitation models of 
both maser and quasi-thermal methanol. It will be partic- 
ularly important to reassess the models of class I methanol 
masers in the light of the new data, since collisions are 
thought to provide the source of energy for excitation. Nev- 
ertheless, for both the class I and class II masers, it is the 
radiative transitions which are decisive for generating the 
population inversions (Cragg et al. 1992). Thus it is not 
surprising that we find the list of transitions which become 
class II masers to be hardly altered by the new collisions. 
What has changed is the details of conditions under which 
individual masers are thermally quenched, and in particu- 
lar the upper limit on gas density in the 6.6-GHz maser 
regime, which has increased by approximately a factor of 10 
to 10 9 cm" 3 . As before, the observational limits suggest the 
majority of 6.6-GHz sources (those without accompanying 
masers at 23.1 or 107.0 GHz) are either near the high den- 
sity limit of maser range (now 10 7 — 10 9 cm -3 ), and/or near 
the low dust temperature threshold (100 < 7d < 150 K). So 
our conclusions regarding the majority of class II methanol 
maser sources are qualitatively unchanged. 

The updated list of maser candidates (Tabled includes 
108 transitions, of which 83 are in common with the corre- 
sponding table of SCG97b. The 14 transitions in the list of 
SCG97b which are absent from Tableware mostly close to 
the threshold for inclusion (7b = 10 4 K). Similarly, a few 
of the additional candidates in this new list are rather close 
to this threshold. The remainder of the new candidates are 
included by virtue of their appearance under the conditions 
of Model 8, suggestive of an intermediate regime spanning 
the switch between class I and class II, which was not exam- 
ined in SCG97b. While the collision model may be further 
improved in the future if rate coefficients for higher J values 
and for ortho-Hb become available, it is encouraging that 
the potentially observable maser lines can be identified irre- 
spective of the knowledge of accurate collisional data. 

It is readily apparent from Table |2] that the warm 
dust/cool gas combination (Models 3 and 7) has the greatest 
potential for producing class II masers. Since this combina- 
tion also produces very intense masers at 6.6 and 12.1 GHz, 
the original modelling of Sobolev & Deguchi (1994) focussed 
on this regime. However, subsequent surveys at several can- 
didate frequencies have shown the weaker maser transitions 
to be rarely detectable. There are 25 sources with masers 
detected in the 3i — 4o A + transition at 107.0 GHz out of 
more than 175 class II methanol maser sources surveyed 
(Val'tts et al. 1995, 1999; Caswell et al. 2000; Minier & 
Booth 2002), as well as a handful of exceptional sources 
with masers known also at several other frequencies. While 
the conditions in these sources are clearly not representa- 
tive of the majority, they provide the opportunity for multi- 
transition fitting to constrain those conditions, as illustrated 
earlier for W3(OH). Here the availability of interferometric 



data makes the results more reliable. Without this we can- 
not be confident that the masers at different frequencies are 
coincident, as is assumed by the modelling. 

The models of the northern maser component in 
W3(OH) are consistent with a warm dust/warm gas regime, 
confirmed above using the new colliisons. The 108.8-GHz 
Oo — l-i E transition is significantly affected by the im- 
proved PFD collision model, not surprisingly since it comes 
from levels so low in energy. We have investigated trends in 
models with parameters in the same vicinity as those de- 
rived previously, without fully accounting for the observed 
upper limit on this maser. Further exploration of the models 
may yield a better fit, but some of these require new phys- 
ical assumptions which are beyond the scope of this paper 
(e.g. different dust models as in Ostrovskii & Sobolev 2002). 
Nevertheless, the detection of strong masers at 19.9, 23.1, 
37.7 and 85.5 GHz in a few sources is evidence of rather 
different conditions elsewhere, for example in NGC 6334F, 
G345.01+1.79, and the southern source in W3(OH). With 
the new collision data there may be new possibilities to ex- 
plain these observations, and we will consider this in more 
detail in future work together with additional observational 
data providing better constraints. 



5 CONCLUSION 

Since the reliability of the model calculations for class II 
methanol masers is influenced by the completeness and ac- 
curacy of the fundamental molecular data employed, we have 
undertaken new model calculations to ascertain firstly the 
effects of including more energy levels, and secondly the 
effects of the new collision rate coefficients. The improved 
molecular data does not lead to dramatic changes in the 
models, as the maser pumping is dominated by radiative 
effects, and the previous models included the essential pro- 
cesses. The list of approximately 100 transitions which are 
potential class II methanol maser candidates is litle changed 
and may be regarded as firmly established. Energy levels up 
to v t —3 and J=18 are found to be sufficient for models with 
dust temperatures up to 300 K. The new collisional rate co- 
efficients of PFD do affect the absolute and relative bright- 
nesses of the masers over the intermediate density range 
10 4 — 10 9 cm -3 . But at gas temperatures above 100 K the 
changes are broadly equivalent to a reduction of the colli- 
sion cross-section, which may be compensated by an increase 
in the gas density. Thus previous conclusions regarding the 
nature of the most common 6.6-GHz class II maser sources 
should be modified to accommodate a slightly higher range 
of gas densities (maximum 10 9 cm -3 ). Fits of the model to 
observations in special sources exhibiting maser emission at 
many frequencies are modified in detail by the new data, but 
the example investigated shows that a comparable fit can be 
obtained under nearby conditions to those previously identi- 
fied. The newly available data therefore allows the modelling 
to proceed with increased confidence. 
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APPENDIX A: METHANOL ENERGY LEVELS 
AND TRANSITIONS 

Fig.Qdisplays energies of E-species methanol up to approx- 
imately 600 K in the ground torsional state «t=0, together 
with the corresponding levels for t; t =l,2,3,4, and for the 
Vt=0 mode of the CO-stretch vibration wco=L In each tor- 
sional or vibrational state 284 rotational levels are included, 
up to an energy threshold corresponding to maximum J=22 
and maximum K=9. The same number of levels are required 
for the A species. The levels are arranged as ladders of in- 
creasing J for each value of K. It is evident that the ladders 
are spaced in an irregular fashion in the torsionally excited 
states, with considerable overlap among the different tor- 
sional states. For example, there is no clear demarcation 
between v t — l and vt=2. The disorder is a result of the large 
torsional-rotational interaction in methanol, which also pro- 
duces irregular spacings in the torsional ground state levels 
(not apparent on this scale). In contrast, the vco=l lev- 
els more closely resemble the vt=0 levels. The i»t=0,l,2,3,4 
energies come from Mekhtiev, Godfrey & Hougen (1999), 
calculated using the Hamiltonian of Xu & Hougen (1995), 
which gives microwave accuracy for the rotational states of 
i>t=0,l,2. The energies of the «co=l levels are calculated 
from the data of Moruzzi et al. (1989), using Taylor se- 
ries fits to observed infrared frequencies, relative to the vi- 
brational ground state energies from Moruzzi et al. (1992). 
The t>co=l levels overlap significantly in energy with both 
v t =3 and v t —4. The v t =l mode of the «co=l vibration (not 
shown) also overlaps in energy with Dt=4, but data is not 
available for the complete range of K required for our cal- 
culations. The extensive study of Moruzzi et al. (1989) also 
includes data on the CH3 rock levels for a single K value in 
the A symmetry species: these are at comparable energy to 
the CO stretch. We may therefore expect that even the large 
number of energy levels illustrated (1704 for each symmetry 
species) will be insufficient under circumstances where the 
i>co=l levels have a significant effect on the maser pumping, 
because levels of comparable energy are omitted. 

The line strengths have been accurately calculated 
for transitions between rotational levels up to J=22 by 
Mekhtiev et al. (1999) for torsional states i> t =0,l,2. This 
includes both transitions within each torsional state, and 
transitions between the different states. We have used ap- 
proximate methods to evaluate the line strengths for the 
higher excited states. These should be sufficient to reveal 
whether or not the higher states are important contribu- 
tors to the maser pumping, and therefore whether the con- 
siderable effort of accurately calcuating more line strengths 
is warranted. Line strengths for AK=0,±1 radiative transi- 
tions within v t =3, and within v t —4, were approximated by 
the corresponding symmetric rotor line strength or Honl- 
London factor for transitions between states of the same 
J and K (e.g. Zare 1988, Table 6.6). Line strengths for 
AK=±1 transitions between v t =3 or «t=4 and the other 
torsional states were set to the product of the symmetric 
rotor line strength with the square of the torsional overlap. 



The torsional overlap between states of different v% and K 
was obtained from the eigenstates of the torsional Hamil- 
tonian of De Lucia et al. (1989) (incorporating the Wang 
transformation for the A-species to obtain the correct ± 
symmetry states). The Einstein A coefficients which gov- 
ern the radiative transition probabilities are calculated from 
the line strengths, and depend also on the square of the 
relevant dipole moment component and the cube of the fre- 
quency. Both the inter- and intra-torsional transitions are 
governed by the permanent dipole moment components: 
Ha = 0.896 Debye for AK=0 transitions and /it, = 1.412 De- 
bye for AK=±1 transitions (Sastry, Lees & Van der Linde 
1981). For transitions among the «co = l levels we used the 
same line strengths as for the corresponding w t =0 transi- 
tions, and the same dipole moment components. Since the 
CO stretch takes place (mainly) along the a-axis, the major 
transitions between vco=l and vco=0 will be transitions 
to v t —0 with AK=0, and these were again approximated by 
the symmetric rotor line strengths. For the CO-stretch tran- 
sition dipole moment we used the value 0.27 Debye obtained 
experimentally by Henningsen & Petersen (1988). 



APPENDIX B: RATE COEFFICIENTS FOR 
COLLISIONS 

The methanol collision partners are assumed to be 20 per- 
cent He and 80 percent para-H2. Pottage et al. (2002, 2004b) 
tabulate rate coefficients for v t =0 levels up to J=9 at 13 tem- 
peratures between 5 and 200 K, for both collision partners 
and both A and E symmetry species of methanol. Note that 
the state designation used by PFD is not exactly the same 
as used here; for example E0 and El refer to E species v t —0 
and Vt = l levels respectively, with AO and Al defined in a 
corresponding fashion. For the A species the sign of K used 
in the most recent tabulations by PFD matches the ± sym- 
metry label. We used linear interpolation between adjacent 
tabulated values for intermediate temperatures, with scal- 
ing by y/T for extrapolation to higher temperatures. When 
there are near-degenerate levels, as occurs for high J and 
K in the A species, the data presented by PFD includes 
rate coefficients for only one state of the degenerate pair; 
the same rate coefficients were used for the other state. The 
rate coefficient for AJ=0=AK asymmetry doublet transi- 
tions between near-degenerate levels was set equal to the 
value for the J=9, K=4 doublet, which is the highest dou- 
blet for which PFD distinguish between the levels. 

Since the method of PFD generates less accurate en- 
ergies than we employ, the states used in the calculation 
of rate coefficients must be mapped onto the set of levels 
used in the excitation calculations. Some care is required 
when the states appear in different energy sequence in the 
two sets, so that an upward transition in one set becomes a 
downward transition in the other, and vice versa. Rate co- 
efficients R for upward and downward transitions between 
any pair of levels i and j obey the detailed balance relation, 
which ensures that a Boltzmann population distribution at 
temperature T is recovered in the high density limit. This 
may be expressed in terms of the following quantity, which 
is independent of the sense (up or down) of the transition. 

giRi^i expCEj/fasT) = gjR^iexpiEi/kBT) 
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where g represents the (2J+1) statistical weight factor, 
E is the energy, and fee is the Boltzmann constant. This in- 
variant expression was used to rescale both the upward and 
downward rates, after ensuring that the energies in both sets 
were referred to the same ground state, in order to account 
for the different energy values attributed to the states in 
our formulation, and specifically for differences in the en- 
ergy sequence. For the remaining transitions where no data 
is available from PFD the rate coefficients for transitions 
downward in energy were defined as described below, with 
the corresponding upward rate coefficients calculated from 
detailed balance. 

For transitions starting from J>9 in Vt~0, the J=9 rate 
coefficient for the transition from the same initial K to the 
same final state was used where possible, extrapolated ac- 
cording to the 1/| A J | propensity rule. When both initial and 
final J are >9, the corresponding J=9 rate coefficient for the 
same AJ and initial and final K values was used, if necessary 
scaled by the 1/|AJ| propensity rule if A J exceeds what is 
available in the PFD data. Since our calculations include 
levels up to either J=18 or J— 22, there are a large number 
of collisional transitions which use these 'extrapolated' rate 
coefficients, but it is not so critical to have accurate values 
as these higher J levels carry only a small portion of the 
population. 

Transitions within v t — l for He interacting with A- 
species methanol were treated the same as v t =0. Data is 
available for T = 20 K only, and was extrapolated to other 
temperatures using the VT rule. The He data was assumed 
in this case to represent para-H2 also. The corresponding 
E-species data was not used as it is insufficiently accurate 
(D. Flower, private communication). 

Rate coefficients for rotational transitions within all 
other torsionally excited states, and within vco — l, were 
based on the propensity rules. These are similar in form 
to those of Lees & Haque, but based on a detailed examini- 
nation of the PFD data for v t =l A species. At temperature 
T = 20 K the rate coefficient 3 x 10~ n em's -1 was used for 
AJ=1, AK=0 transitions, together with 1 x 10 -11 em's -1 
for AJ=0, |AK|>0 transitions. In both cases 1/|AJ| scaling 
was used for larger values of AJ. The value 3 x 10~ n em's -1 
was used for the A species AJ=0=AK asymmetry doublet 
transitions. Rate coefficients at different values of tempera- 
ture were scaled by \/T. 

Pottage et al. (2004a) calculated cross-sections for tran- 
sitions between v%=\ and «t=0, but not rate coefficients. 
Cross-sections for torsionally inelastic transitions were found 
to be typically two orders of magnitude smaller than cross- 
sections for torsionally elastic transitions. Furthermore they 
exhibited no clear propensity for particular AJ or AK. Based 
on the typical magnitudes of the v t =0 rate coefficients, we 
adopted the value 5 x 10 -14 em's -1 for all downward tran- 
sitions involving changes in v t or vco ■ 

REFERENCES 

Castor J. I., 1970, MNRAS, 149, 111 

Caswell J.L., Yi J., Booth R.S., Cragg D.M., 2000, MNRAS, 313, 
599 

Cragg D.M., Johns K.P., Godfrey P.D., Brown R.D., 1992, MN- 
RAS, 259, 203 



Cragg D.M., Sobolev A.M., Ellingsen S.P., Caswell J.L., Godfrey 

P.D., Salii S.V., Dodson R.G., 2001, MNRAS, 323, 939 
Cragg D.M., Sobolev A.M., Caswell J.L., Ellingsen S.P., Godfrey 

P.D., 2004, MNRAS, 351, 1327 
De Buizer J.M., Pina R.K., Telesco CM., 2000, ApJS, 130, 437 
De Lucia F.C., Herbst E., Anderson T., Helminger P., 1989, 

J.Mol.Spec, 134, 395 
Ellingsen S.P., von Bibra M.L., McCulloch P.M., Norris R.P., 

Deshpande A.A., Phillips C.J., 1996, MNRAS, 280, 378 
Ellingsen S.P., Cragg D.M., Lovell J.E.J., Sobolev A.M., Rams- 
dale P.D., Godfrey P.D., 2004, MNRAS, 354, 401 
Haschick A.D., Baan W.A., Menten K.M., 1989, ApJ, 346, 330 
Henningsen J.O., Petersen J.C., 1988, J.Opt.Soc.Am. B, 5, 1848 
Lees R.M., Baker J.G., 1968, J.Chem.Phys., 48, 5299 
Lees R.M., Haque S.S., 1974, Can.J.Phys., 52, 2250 
Lovas F.J., 2004, J.Phys.Chem.Ref.Data, 33, 177 
Malyshev A.V., Sobolev A.M., 2003, Astro. Astrophys Tr., 22, 2 
Mekhtiev M.A., Godfrey P.D., Hougen J.T., 1999, J.Mol.Spec, 
194, 171 

Menten K.M., Johnston K.J., Wadiak E.J., Walmsley CM., Wil- 
son T.L., 1988, ApJ, 331, L41 

Menten K.M., Reid M.J., Pratap P., Moran J.M., Wilson T.L., 
1992, ApJ, 401, L39 

Minier V., Booth R.S., 2002, A&A, 387, 179 

Minier V., Booth R.S., Conway J.E., 2000, A&A, 362, 1093 

Moruzzi C, Strumia F., Carncsecchi P., Lees R.M., Mukhopad- 
hyay I., Johns J.W.C, 1989, Infrared Phys., 29, 583 

Moruzzi C, Strumia F., Moracs J.C.S., Lees R.M., Mukhopad- 
hyay I., Johns J.W.C, Winnewisser B.P., Winnewisser M., 
1992, J.Mol.Spec., 153, 511 

Moscadclli L., Menten K.M., Walmsley CM., Reid M.J., 1999, 
ApJ, 519, 244 

Norris R.P., Whiteoak J.B., Caswell J.L., Wieringa M.H., Cough 

R.G., 1993, ApJ, 412, 222 
Ostrovskii A.B., Sobolev A.M., 2002, in Migenes V., Reid M.J., 

eds, Proc. IAU Symp. 206, Cosmic Masers: from Protostars 

to Black Holes., Astron.Soc.Pac, San Francisco, p. 183 
Peng R.S., Whiteoak J.B., 1993, MNRAS, 260, 529 
Pestalozzi M., Minier V., Booth R.S., 2004, A&A, in press 
Pottage J.T., Flower D.R., Davis S.L., 2001, J Phys. B, 34, 3313 
Pottage J.T., Flower D.R., Davis S.L., 2002, J Phys. B, 35, 2541 
Pottage J.T., Flower D.R., Davis S.L., 2004a, J Phys. B, 37, 165 
Pottage J.T., Flower D.R., Davis S.L., 2004b, MNRAS, 352, 39 
Sastry K. V.L.N. , Lees R.M., Van der Linde J., 1981, J.Mol.Spec, 

99, 228 

Slysh V.I., Kalcnskii S.V., Val'tts I.E., 1993, ApJ, 413, L133 
Slysh V.I., Kalcnskii S.V., Val'tts I.E., 1995, ApJ, 442, 668 
Slysh V.I., Val'tts I.E., Kalcnskii S.V., Larionov CM., 1999, Astr. 

Reports, 43, 657 
Sobolev A.M., Deguchi S., 1994, A&A, 291, 569 
Sobolev A.M., Cragg D.M., Godfrey P.D., 1997a, A&A, 324, 211 
Sobolev A.M., Cragg D.M., Godfrey P.D., 1997b, MNRAS, 288, 

L39 (SCG97b) 

Sutton E.C, Sobolev A.M., Ellingsen S.P., Cragg D.M., 
Mehringer D.M., Ostrovskii A.B., Godfrey P.D., 2001, ApJ, 
554, 173 

Val'tts I.E., Dzura A.M., Kalcnskii S.V., Slysh V.I., Booth R.S., 

Winnberg A., 1995, A&A, 294, 825 
Val'tts I.E., Ellingsen S.P., Slysh V.I., Kalcnskii S.V., Otrupcck 

R., Voronkov M.A., 1999, MNRAS, 310, 1077 
Wilson T.L., Hiittcmcistcr S., Dahmen C, Henkel C, 1993, A&A, 

268, 249 

Wilson T.L., Walmsley CM., Menten K.M., Hermscn W., 1985, 
A&A, 147, L19 

Voronkov M.A., Sobolev A.M., Ellingsen S.P., Ostrovskii A.B., 

Alakoz A., 2004, Astrophys. Space Sci., in press 
Xu L.H., Hougen J.T., 1995, J.Mol.Spec, 173, 540 
Zare R.N., 1988, 'Angular Momentum', Wiley, New York 



